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Abstract 

The discontinuity of the vorticity is written as a function of the vector Tgrads, 
(where T is the temperature and s the specific entropy). The expression is obtained 
thanks to potential equations and independently of the mass conservation and the 
equation of momentum balance. 
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1 Introduction 



The aim of this note is to prove that, in the most general instationary case of 
perfect compressible fluids, across a shock wave we have the relations: 

[w(rot v) t g] = n A [T grad s] 
[(rot v) n ] = 

where v is the velocity vector of the fluid, u is the fluid velocity with respect 
to the shock wave, T is the temperature, s is the specific entropy, the indices 
tg and n indicate the tangential and normal components to the shock wave of 
the vector rot v and the discontinuity of a tensorial quantity a is denoted by 
[a] (see[l]). 
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We take into account the following shock conditions: 

[v tfl ] = 

■1 



u -\- h 







(2) 
(3) 



dp 

where h denotes the specific enthalpy (dh = Tds H ). 

P 

We do not use the following shock conditions: 

H = o 



p + pu z 







(4) 
(5) 



where p is the density and p the pression of the fluid. 



We will use the potential equations obtained by P. Casal [2] or J. Serrin [3] 
and expressing another form of the equations of the motions of compressible 
perfect fluids. 

The motion of a compressible fluid is represented by a continuous mapping of 
a reference three-dimensional space D Q in the physical space D t occupied by 
the fluid at time t: 

x = (p t (X), x G D t , X G D , 

or equivalently by a continuous mapping $ : W D — > W, z = $(Z) where W 
is a four-dimensional reference space and W the physical time-space, 



G W and Z 




G W a . 



We assume the motion has a shock wave localized on a surface S(t) propagating 
in D t , image by tp t of a surface S Q (t) propagating in D Q . We denote by n Q and 
n the unit normal vectors to S Q (t) and S(t) respectively, and g Q and g their 
respective velocities; then u = n T v — g, where T denotes the transposition. 

Equivalently, S Q and S are respectively the two corresponding surfaces prop- 
agating in W and W; N G and N are the associated normal vectors, 

N T =(-g,n T ). 

Consequently, $ is a differential mapping on W , except on S Q ; its Jacobian 
matrix is denoted by dz/dZ, F denotes the Jacobian matrix of ip t : 



rfx = vdt + FrfX. 
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2 Exterior derivative 



The covector C = v T has an inverse image C Q in D a such that 



C = CF 



The exterior derivative of the form C is a 2-form which is isomorph to the 
vector rot v. It the image of the 2-form which is the exterior derivative of C Q 
isomorph to the vector rot Q [C T ): 



rot G is the rotational on D Q (see reference [4]). 

The discontinuity of the vorticity vector comes from to parts: one part comes 
from the discontinuity of its image rot Q (C T ) , and the other part comes from 
the discontinuity of the Jacobian F. 



3 Discontinuity of the Jacobian F 

dz/dZ is a linear mapping transforming any tangent vector to S a (t) in a 
tangent vector to S(t). If we denote by n Q ' = —n o /g 0: we obtain: 



where indices 1,2 indicate quantities upstream and downstream the shock. 
Consequently, 





(6) 



[F\ = [vj n G , n Q =n 1 — = n 1 — 



n 







Taking into account Eq. (2) we obtain 



[v] = [u] n, 



(7) 



[F] = [u]nn ' T , 



(8) 





where I is the identity matrix. 
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4 Discontinuity of the 1-form C Q 



Due to the fact that C Q = v T F and by using Eq. (7) and Eq. (8), we obtain: 

[C ]=[u 2 + gu]n ' T . (10) 

5 Discontinuity of potentials 



The dot denotes the material derivative, Q is the body force potential. We 
consider the two quantities <p(t, X) and ip(t, X) (denoted potentials) such that 

<p = 0(t,X), (11) 

^ = 7 (f,X). (12) 

(3 and 7 are two scalar fields defined in each point of the flow and such that 

P(t,X) = ^v 2 -h-n, 7 (t,X) = T. 

There exists a covector B function only of X such that 

We can verify that Eq. (13) together with Eq. (11) and Eq. (12) are equivalent 
to potential equations proposed by P. Casal in [2] and J. Serrin in [3]. 
With the condition of adiabaticity 

s = 0, 

and the equation of balance of mass 

^ + div(pv) = 0, 
we obtain the complete set of motion equations. 

We can choose ip and ip null on the shock wave and continuous through the 
shock surface in the following manner: 

*=L p{T ' x)dT ' ^ = U {T ' x)dT < (i4) 

where t = /(X) is the equation of the shock surface S a (t) which is assumed 
regular. 
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Potential ip being continuous through the shock surface, 



dp 
dX 



[<f\ no'' 



-w 2 -h 



n r 



By using Eq. (2), Eq. (3) and Eq. (10) we get: 

dp 



dX 



[Co] 



(15) 



Due to the fact that ip is null on the shock wave, Eq. (13) expresses B is 
continuous through the shock: 



[B] = 



(16) 



6 Discontinuity of the image of the vorticity 



Let us consider 



dp 



W = Cj - ^ J , then rot W G = rot G C G T . 

Eq. (13) yields W Q = Tgrad Q s, the value of grad Q s being defined on W a . 
Due to Eq. (15), W Q is continuous through the shock and we get: 
'dW 



dX 



n Q ' , [rot W ] = n Q ' A 



rot C T ] = n>[T grade s]. (17) 

The discontinuity of the image of the vorticity is only tangential. Using the 
previous results and application <3>, we can verify that this property of the 
vorticity remains true. 



7 Discontinuity of the vorticity 



The results obtained by Eqs (6), (9) and (17) allow to obtain Formulae (1). 
This expression general for non stationary perfect compressible fluids is dif- 
ferent from the result given by Hayes [5]. This is due to the fact the result is 
obtained thanks to Eq. (3) of conservation of energy. It neither uses Eq. (5) of 
the balance of the quantity of motion nor Eq. (4) of the conservation of mass. 
To obtain Formulae (1), the knowledge of the enthalpy field is only necessary. 
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In the special case of a stationary, iso-energetic, irrotational motion upstream 
of the shock, the relation can be expressed with the help of the curvature 
tensor of the shock surface [6]. 
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